Dynamics and symmetries of flow reversals in turbulent convection 
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Based on direct numerical simulations and symmetry arguments, we show that the large-scale 
Fourier modes are useful tools to describe the flow structures and dynamics of flow reversals in 
Rayleigh-Benard convection (RBC). We observe that during the reversals, the amplitude of one of 
the large-scale modes vanishes, while another mode rises sharply, very similar to the "cessation-led" 
reversals observed earlier in experiments and numerical simulations. We find anomalous fluctuations 
in the Nusselt number during the reversals. Using the structures of the RBC equations in the Fourier 
space, we deduce two symmetry transformations that leave the equations invariant. These symmetry 
transformations help us in identifying the reversing and non-reversing Fourier modes. 
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Many experiments [TH7] and numerical simulations [SJ 
I8HTU] on turbulent convection reveal that the velocity 
field of the system reverses randomly in time (also see re- 
view articles [H]). This phenomenon, known as "flow re- 
versal" , remains ill understood. This process gains prac- 
tical importance due to its similarities with the magnetic 
field reversals in geodynamo and solar dynamo [12]. In 
this letter, we study the dynamics and symmetries of flow 
reversals in turbulent convection using the large-scale 
Fourier modes of the velocity and temperature fields. 

The experiments and simulations performed to explore 
the nature of flow reversals are typically for an ideal- 
ized convective system called Rayleigh-Benard convec- 
tion (RBC) in which a fluid confined between two plates 
is heated from below and cooled at the top. Detailed 
measurements show that the first Fourier mode vanishes 
abruptly during some reversals [3j 2]. These reversals 
are referred to as "cessation- led" . Recently Sugiyama et 
al. [5] performed RBC experiments on water in a quasi 
two-dimensional box, and observed flow reversals with 
the flow profile dominated by a diagonal large-scale roll 
and two smaller secondary rolls at the corners. They at- 
tribute the flow reversals to the growth of the two smaller 
corner rolls as a result of plume detachments from the 
boundary layers. 

Several theoretical studies performed to understand 
reversals in RBC provide important clues. Broadly, 
these works involve either stochasticity (e.g., "stochas- 
tic resonance" [8j [13]), or low-dimensional models with 
noise [I4j [15] . Mishra et al. [10] studied the large-scale 
modes of RBC in a cylindrical geometry and showed 
that the dipolar mode decreases in amplitude and the 
quadrupolar mode increases during the cessation-led re- 
versals. Regarding dynamo, low-dimensional models con- 
structed using the large-scale modes and symmetry argu- 
ments reproduced dynamo reversals successfully [TJ [16] . 

The theoretical models described above only focus on 
the large-scale modes. Here too, they provide limited 
information about these modes due to small number of 
measuring probes. In this letter we compute large-scale 



and intermediate-scale Fourier modes accurately using 
the complete flow profile. This helps us in quantita- 
tive understanding of the dynamics and symmetries of 
the RBC system. We also show that these modes can 
describe the diagonal and corner rolls of Sugiyama et 
al. [5] , as well as the cessation- led reversals of Brown and 
Ahlers [3 . The properties of the modes for convection 
are contrasted with those for dynamo. 

The equations governing two-dimensional Rayleigh- 
Benard convection under Boussinesq approximation are 

+ (u • V)u = -VP + PrV 2 u + RaPrTy, (1) 

ot 

^ + (u-V)T = V 2 T, (2) 

V • u = 0, (3) 

where u is the velocity field, T is the temperature field, P 
is the pressure, and y is the buoyancy direction. The two 
nondimensional parameters are the Prandtl number Pr, 
the ratio of the kinematic viscosity v and the thermal dif- 
fusivity ft, and the Rayleigh number Ra = agAd 3 /(vK), 
where a is the thermal expansion coefficient, d is the 
distance between the two plates, A is the temperature 
difference between the plates, and g is the acceleration 
due to gravity. The above equations have been nondi- 
mensionalized using d as the length scale, the thermal 
diffusive time d 2 / n as the time scale, and A as the tem- 
perature scale. 

We solve the above equations in a closed box geometry 
of aspect ratio Y = 1 (denoted by Pi) and T = 2 (de- 
noted by P2) using Nek5000 [17] . an open source spectral- 
element code. We apply no-slip boundary condition on 
all the walls. The top and bottom walls are assumed 
to be perfectly conducting, while the side walls are as- 
sumed to be insulating. We use 28 x 28 spectral elements 
for Pi and 48 x 28 spectral elements for P2, and 7 th or- 
der polynomials for resolution inside the elements. Thus, 
the effective grid resolution for the Pi and P2 boxes are 
196 x 196 and 336 x 196 respectively. The concentration 
of grid points is higher near the boundaries in order to 




FIG. 1. Top panel: Convective flow profile computed from a spectral element simulation for a box of aspect ratio r = 1, Prandtl 
number Pr = 1. Rayleigh number Ra = 2 x 10 7 : (a) before the reversal (t = 10.52 thermal diffusive time units), (b) during 
the reversal (t = 10.68), (c) after the reversal (t = 10.79). Bottom panel: Velocity field corresponding to (a) ui,i, (b)u2,2, (c) 
— ui,i. The flow patterns of the bottom panel are good approximations of the corresponding patterns of the top panel. The 
corner rolls of (a,c) are obtained as a result of superposition of (d,f) with (e). 



resolve the boundary layer. We perform our simulations 
for Ra = 2 x 10 7 , 10 8 , 10 9 for B u and for Ra = 10 7 , 
2 x 10 7 , 10 8 for £?2, till several thermal diffusive time 
units. Pr = 1 for all our runs. We observe flow reversals 
only for Ra = 2 x 10 7 (#i) and Ra = 10 7 (B 2 ). These 
results are consistent with those of Sugiyama et al. [5] . 

In Fig. [TJa,b,c) we display three frames of the velocity 
and temperature fields for the Ra = 2 x 10 7 run for the 
T = 1 box. The three frames illustrate the flow profiles 
before the reversal ((a), at t = 10.52 thermal diffusive 
time units), during the reversal ((b), at t = 10.68), and 
after the reversal ((c), at t = 10.79) (see videos in [18). 
They are very similar to those presented by Sugiyama et 
al. 0. To gain further insights into the reversal dynam- 
ics, we decompose the velocity and temperature fields 
into Fourier modes 

u = ^ w m ,n sin(mA: c x) cos(ran/), (4) 

m,n 

v = # m ,n cos(mk c x) sin (ran/), (5) 

m,n 

T = ^2 Tm,n cos(mk c x) sin(rar?/) , (6) 



where u = (u,v) is the velocity field, k c = 7r for T = 1, 
and k c = tt/2 for T = 2. 

We compute the Fourier amplitudes using FFTW li- 
brary through PyFFTW interface [19] . For these trans- 
forms we interpolate the Nek5000 data to a uniform 
128 x 128 grid. In Fig. []Jd,e) we display the velocity 
profiles of the primary modes k = (1, 1) and (2, 2), which 
are a single roll, and four rolls respectively. The sign of 
the Fourier mode of Fig. [TJf) is reversed compared to 
that of Fig. [ljd) . The diagonally oriented roll and the 
corner rolls of Fig. [lja,c) can be well approximated as a 
superposition of Fig. [TJd,f) and Fig. [TJe) with appropri- 
ate amplitudes. The corresponding primary modes for 
T = 2 box are k = (2, 1) and (2, 2). 

In Fig. [5] we plot time series of (a) the vertical veloc- 
ity V y at (0.25,0.25), (b) the dominant Fourier modes 
^1,1,^2,2, (c) the ratio |^2, 2/^1,1 1, and (d) the Nusselt 
number, for T = 1. The time series indicates that the 
mode k = (1, 1) dominates the mode (2,2) between the 
reversals. During the reversals however, the mode v\ t \ 
crosses zero, while £2,2 shows a spike. The mode 
overshoots around 40% before it attains the steady-state. 
Vanishing of £>i 5 i, and the peaking of £2,2 (the corner 
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FIG. 2. For T = 1, Pr = 1, and Ra = 2xl0 7 , time series of (a) 
y component of the real velocity field at the point (0.25, 0.25), 
(b) the Fourier amplitudes £1,1 (blue) and 7)2,2 (red), (c) the 
ratio I £2, 2/^1,1 1, and (d) Nusselt number (Nu). During the 
reversal, the sign of mode #1,1 changes sign after vanishing, 
while the mode 62,2 rises sharply. The Nusselt number shows 
large fluctuations during the reversals. 



rolls) are the reasons for the four rolls appearing during 
the reversal (Fig. [TJb)). This phenomenon is the same as 
the cessation-led flow reversals reported by Brown and 
Ahlers [3] and Mishra et al. [TO] . 

The mode (1,1) changes sign due to the reversals but 
the mode (2, 2) does not. This is due to certain sym- 
metries obeyed by the governing equations which will be 
discussed later. As a result, the velocity field near the 
corners does not reverse (see Fig. [T]). Fig. [2] also shows 
that the Nusselt number becomes negative during the re- 
versals, i.e., the fluid loses heat energy to the plates for 
a small time interval while the flow reverses. 

For the Y = 2 box, the dynamics of the modes are ex- 
actly the same as above, except that the Fourier mode 
(2, 1) takes the role of the (1, 1) mode. That is, the pri- 
mary modes for the Y = 2 box are (2, 1) and (2, 2), which 
represent the two large rolls, and the corner rolls respec- 
tively of Fig. [3] Note however that in the simulations 
of Brueur and Hansen [9 for Y — 2 and Pr — » 00 fluid, 
the most dominant modes were (1,1) and (2,1). Thus, 
a variation of the Prandtl number for the same box can 
change the flow pattern. For our simulations the aver- 
age intervals between consecutive reversals for Y = 1 and 
T = 2 are approximately 0.6 and 0.05 thermal diffusive 
time units respectively. The duration of the reversals is 
approximately 0.03 time units for both the cases. 

The flow reversals are not observed for Y = 1 and Ra = 
10 8 , 10 9 , and for Y = 2 and Ra = 2 x 10 7 , 10 8 till several 
thermal diffusive time units (maximum 10), though they 
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FIG. 3. Convective flow profile for Y — 2, Pr — 1, (a) re- 
versing case Ra — 10 7 , and (b) non-reversing case Ra — 10 8 . 
The corner rolls of subfigure (a) are much larger than those of 
subfigure (b). The ratios of the time-averaged Fourier ampli- 
tudes (\i)2,2 \}t/(\v2,i \)t for the cases are 0.12(a) and 0.04(b). 
Small amplitude of £2,2 makes the growth of the corner rolls 
difficult. 



may occur after some more time, a result consistent with 
that of Sugiyama et al [5]. As the Rayleigh number is 
increased, the flow reversals become more difficult due 
to suppression of the corner rolls ((2,2) mode) by the 
dominant roll structure, which is quantified by the ratio 
of the time- averaged Fourier amplitudes of these modes. 
For our simulations on Y = 2, the ratio (1^2,2 |)t/(|^2,i|)t 
ranges from 0.12 for Ra = 10 7 (reversing) to 0.04 for 
Ra = 10 8 (non-reversing), which is consistent with the 
convective flow profiles exhibited in Fig. |3|a,b) for these 
two cases. For Y = 1, the corresponding ratio ranges 
from 0.44 for Ra = 2 x 10 7 to 0.12 for Ra = 10 9 . 

We deduce interesting features on the generation and 
symmetry of the Fourier modes using the structure of the 
two-dimensional RBC equations in the Fourier space [20] : 



dui(k) 
dt 



flfl(k) 
dt 



~ ik o ^2 Uj(p)ui(q) - ikiP(k) 

k=p+q 

+RaPr§(k)S it2 - Prk 2 Ui(k) 



(7) 



-ikj %(p)#(q)+^(k)-/c 2 #(k)(8) 

k=p+q 



where is the perturbation of the temperature about 
the conduction state. For a two-dimensional box, the 
Fourier modes (coefficients of the basis functions) of the 
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fields are of the type E = (even, even), = (odd, odd), 
M eo = (even, odd), and M oe = (odd, even), where we 
refer to E : and M as even, odd and mixed modes 
respectively. The nonlinear interactions generate new 
modes with k = p + q. For example, modes (rai,ni) 
and (772,2,^2) generate modes (|rai ±7712!, |ni ±712 1) (note 
that sin(rax) is a superposition of both k x = ±m). The 
nonlinear interactions satisfy the following properties: 
O x O = E; O x E = O; E x E = E; E x M e0y0e = 
M e0)0e ; O x M e0j0e = M oe , eo ; M e0j0e x M e0)0e = E\ and 
M e0i0e x M oe ^ eo = O. Here O x O = E means that two 
O-modes interact to yield an E-mode. As a result of 
these properties, the two symmetry operations that keep 
Eqs. ([7| [8| invariant are: 

1. (E -> £, O -O, M eo M eo , M oe -> -M oc ) 

2. (E -> O -O, M eo -M eo , M oe -> M oe ) 

That is, for case (i), if {E, O, M eo , M oe } is a solution of 
Eqs. ([7||8|, then {£", — O, M eo , — M oe } is also a solution of 
these equations. The system explores the solution space 
allowed by these symmetry properties. Note that the 
above properties are universal, that is, they are indepen- 
dent of the box geometry, Prandtl number, etc. 

To relate our simulations with the above mentioned 
symmetry, we observe that the (1,1) (of O-type) and 
(2, 2) (of E-type) modes are the primary modes for the 
r = 1 system. Our simulations show that the mode 
(1,1) and other O-type modes switch sign after a rever- 
sal, while the sign of (2, 2) and other i^-type modes re- 
mains unchanged. The M-type modes have very small 
energy, hence both the symmetries reduce to one, which 
is the symmetry of the T = 1 system. For the T = 2 box, 
the primary modes are (2, 1) (of M eo -type) and (2, 2) (of 
i^-type). These primary modes and subsequently gen- 
erated modes satisfy the second symmetry mentioned 
above since the M eo -type modes flip, while the E-type 
modes do not. By analogy, we expect the solution of a 
T = 1/2 box (with finite Pr) to satisfy the first symmetry, 
while the modes of Breuer and Hansen [9] satisfy the sec- 
ond symmetry (with (1,1) and (2, 1) as primary modes). 
These symmetry arguments are general, and thus useful 
for understanding the dynamics of the Fourier modes. 

We also point out that in magnetohydrodynamics 
(MHD), {u u,b -b} is a symmetry of the MHD 
equations, so all the Fourier modes of the magnetic field 
b change sign after the reversal. The dynamical equa- 
tions of RBC do not have such global symmetry. This 
is one of the critical differences between flow reversals of 
RBC and magnetic field reversals of dynamo. 

To conclude, our numerical simulations and symmetry 
arguments of RBC demonstrate the usefulness of large- 
scale Fourier modes in describing various features of flow 
reversals, such as the large-scale diagonal roll and corner 
rolls of Sugiyama et al [5], and the cessation- led reversals 



reported by Brown and Ahlers [3] and Mishra et al. [TO] . 
We also find that the Nusselt number fluctuates wildly 
during the flow reversals. We exploit the structures of the 
RBC equations in the Fourier space to identify its sym- 
metry transformations. The symmetry arguments and 
the reversal dynamics described in this paper are quite 
general, and they would be useful in understanding of 
reversals in three-dimensional convective flow as well as 
in dynamo. 
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